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to  the  farfield  ranges,  as  long  as  the  square  of  the  index  of  refraction 
is  large  and  the  Sommerfeld  numerical  distance  is  small.  For  the  sub- 
surface-to-air,  air-to-subsurface ,  and  subsurface-to-subsurface  propagation 
cases,  the  additional  restriction  that  the  measurement  distance  be  greater 
than  three  times  the  burial  depth  of  the  source  and/or  receiver  must  be 
met. 
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IMAGE-THEORY  ELECTROMAGNETIC  FIELDS  OF  HORIZONTAL  DIPOLE 
ANTENNAS  IN  PRESENCE  OF  CONDUCTING  HALF-SPACE 

INTRODUCTION 


During  the  past  several  years,  finitely  conducting  earth-image  theory 
techniques  have  proved  quite  useful  in  determining  the  quasi-static  fields  of 
antennas  located  near  the  earth's  surface  for  both  single- layered  and  multi¬ 
layered  earths.  (For  detailed  references,  see  Bannister . * »2)  The  quasi-static 
range  is  defined  as  that  range  where  the  measurement  distance  is  much  less  than 
a  free-space  wavelength. 

Physically,  the  essence  of  the  quasi-static  range  finitely  conducting 
earth-image  theory  technique  is  to  replace  the  finitely  conducting  earth  by  a 
perfectly  conducting  earth  located  at  the  (complex)  depth  d/2,  where  d  =  2/y: 
and  y i  =  [ioiUgCoj  +  iaiEj)]1/2  is  the  propagation  constant  in  the  earth.  (See 
figure  1  for  the  image-theory  geometry.)  Analytically,  this  corresponds  to 
replacing  the  algebraic  "reflection  coefficient,"  (Uj  -  X)/(u1  +  X) ,  in  the 
exact  integral  expressions  by  exp(-Xd),  where  X  is  the  variable  of  integra¬ 
tion.  3  For  antennas  located  at  or  above  the  earth's  surface,  the  general 
image-theory  approximation  is  valid  throughout  the  quasi-static  range. 1.2 


IMAGE 


Figure  1.  Image-Theory  Geometry 
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The  major  disadvantage  of  the  finitely  conducting  earth-image  theory 
technique  is  that  mainly  it  has  been  applied  only  in  the  quasi -static  range. 
However,  because  this  range  includes  the  critical  launching  area,  one  would 
expect  to  be  able  to  extend  it  into  the  intermediate  range  where  the  principal 
field  propagation  proceeds  as  though  it  were  over  a  perfectly  conducting  plane. 

Recently^.4  we  have  shown,  for  horizontally  polarized  sources,  that 
finitely  conducting  earth-image  theory  techniques  are  not  limited  to  the  qua¬ 
si-static  range  alone.  That  is,  by  replacing  the  horizontally  polarized 
algebraic  "reflection  coefficient,"  (Uj  -  Uq)/^  +  uQ)  ,  by  exp(-uQd)  ,  we 
demonstrated  that  finitely  conducting  earth-image  theory  techniques  can  be 
utilized  at  any  range  from  the  source.  Mohsen^  has  validated  and  extended 
these  results  to  include  higher  order  terms  that  correspond  to  multiple  images 
at  the  same  location.  Mahmoud  and  Metwally,^  employing  discrete  and  discrete- 
plus-continuous  images,  have  computed  satisfactorily  the  change  in  the  input 
impedance  of  a  vertical  magnetic  dipole  (VMD)  due  to  the  presence  of  the  earth. 

It  is  the  purpose  of  this  report  to  show  that  nearfield  and  farfield 
range  finitely  conducting  earth-image  theory  techniques  also  can  be  employed 
for  determining  the  fields  produced  by  horizontal  electric  dipole  (HED)  and 
horizontal  magnetic  dipole  (HMD)  antennas  (which  are  a  combination  of  verti¬ 
cally  and  horizontally  polarized  sources). 

In  this  report,  the  HED  and  HMD  are  situated  at  height  h  with  respect  to 
a  cylindrical  coordinate  system  (p,<fi,z)  and  are  assumed  to  carry  a  constant 
current  I.  The  HED  (of  infinitesimal  length  t)  is  oriented  in  the  x  direction 
while  the  axis  of  the  HMD  (of  infinitesimal  area  A)  is  oriented  in  the  y 
direction.  The  earth,  which  is  assumed  to  be  a  homogeneous  medium  with  con¬ 
ductivity  Oj  and  dielectric  constant  £j(=  ere0) ,  occupies  the  lower  half-space 
(z  <  0)  and  the  air  occupies  the  upper  half-space  (z  >  0) .  The  magnetic  per¬ 
meability  of  the  earth  is  assumed  to  equal  u0,  the  permeability  of  free  space. 
Meter-kilogram-second  (MKS)  units  are  employed  and  a  suppressed  time  factor  of 
exp(iwt)  is  assumed. 


LOCATION  OF  IMAGE  DEPTH 


Basic  antenna  theory  tells  us  that  the  fields  produced  by  a  current-car¬ 
rying  wire  of  any  length,  when  placed  over  a  perfectly  conducting  earth,  can 
be  represented  by  the  combined  fields  of  the  wire  and  its  image. 7  if  the 
finitely  conducting  earth  could  be  replaced  with  a  perfectly  conducting  earth 
at  some  specified  depth  below  the  surface  of  the  finitely  conducting  earth, 
we  then  could  use  standard  image  theory  to  locate  the  antenna  image  and  the 
resulting  fields. 

Several  methods  are  available  for  deriving  the  depth  of  a  perfectly  con¬ 
ducting  plane  that  can  be  used  to  replace  a  finitely  conducting  earth.  One  of 
the  most  general  methods**  is  to  equate  the  wave  impedances  for  grazing  inci¬ 
dence  at  the  surface  (z  =  0)  for  the  two  cases  shown  in  figure  2. 
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CASE  B 


z  =  0 
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Figure  2.  Replacement  of  a  Finitely  Conducting  Earth  With  a 
Perfectly  Conducting  Earth  at  Depth  z ^ 

For  case  A,  for  transverse  electric  (TE)  propagation, 

za  V1  -  '^y\  ’ 


where  Yq  =  -w2u0e0,  y 2  =  iuuigCaj  +  iioe^,  and 


/  i^o  \1/2 
nl  "  laj  +■  iaje lJ 


For  case  B,  we  can  write 


ZB  '  n0  tanhCYot!)  , 


where  n0  =  \ju0/^0  ~  120ir. 


For  small  values  of  zi  (i.e.,  Iy^I  i0.5),  tanh(Y0Zj)  -  y0zx»  and 

ZB  ‘  Vozl  s  i(UWozl  *  (4) 

Equating  the  two  impedances  results  in 

Z1  '  - - /."■1  ?  -  ~r==  •  ,  ■  C5) 


Y1  V1  ‘  Y0/Yi  VY 1  '  YC 
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Since  the  image  depth  is  equal  to  2z^t  we  see  that,  for  TE  propagation, 
the  image  depth  d^  for  a  wire  on  the  surface  of  a  finitely  conducting  earth 
can  be  expressed  as 


2hx 


TE  yjyl  "  Y0  Vl  -  l/n2  ' 

subject  to  the  condition  that  | Yq 2 1 1  1.  0.5  (i.e.,  |  yjn2  -  1 1  2,  where  n2  = 

yVyD- 


(6) 


Similarly,  for  transverse  magnetic  (TM)  propagation,  since 

ZA  =  nl  V1  '  yl/y\  ’ 


(7) 


then, 


dTM  ’  =  fW1  *  i/"2  *  C8) 

subject  to  the  condition  that  [ yqz x |  <_  0.5  (i.e.,  |n2/  \nz  -  1|  >_  2) . 

For  normal  incidence,  or  if  |n2|  >  15,  equations  (6)  and  (8)  reduce  to 
the  well-known  result1 ,3 ,8 


^TE  ~  ^TM  ”  d  '  2/y^  , 


(9) 


where  y1 
Oj  <<  we 


'  Viuw oal 
oer‘ 


for  Oj 


>>  we0er  and  Yj 


SOlTOj/  Vs  ia,VlJ0£0er  f°r 


Another  way  to  determine  the  finitely  conducting  earth-image  depth  is  to 
compare  the  results  obtained  from  image  theory  with  known  analytical  results. 
For  z  =  h  =  0,  the  HED  Hertz  vector  is  exactly  equal  to9 


n  = 

X 


U 


4iriwe0  (y\  -  Yg) P " 


(1  +  YnP)e 


"Y  nP 


(1 


-Y  p' 

+  Y,p)e  l 


1  • 


(10) 


The  image-theory  result  is 


u 

[.“V  e'Vi] 

4iriw€g 

.  p  pi  . 

where 


(11) 


When  ReYjP  >>  1,  equation  (10)  reduces  to 


! 


'  * 

;  i 

t 

T - 

A _ 
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IT 


U 


x  4TTiuef,p3 


(f  -  ^  * 


Y0P)e 


-YqP 


while  equation  (11)  becomes 


It 


x  47riue0p3 


U /(1  * 


-YnP 

Y0P)e  0  • 


(12) 


(13) 


Equating  equations  (12)  and  (13)  results  in 


d 


TE 


(14) 


which  is  identical  to  equation  (6) . 

Following  Wait  and  Spies, ^  another  way  to  determine  the  image  depth  is 
to  expand  the  function 


in  a  Taylor  series  about  uQ  =  0,  resulting  in 


where  d  is  given  by  equation  (6) . 

The  introduction  of  exp(-uQd)  into  the  n  integral  equation  yields  an 

image  at  a  distance  h  +  d  from  the  earth's  surface  (see  figure  1),  while 
higher  order  terms  would  correspond  to  multiple  images  at  the  same  location.  * 


HED  AIR-TO-AIR  PROPAGATION 


When  h  and  z  are  >  0,  the  Sommerfeld  integral  expressions  for  the  HED 


Hertz  vector  are 


9-11  “ 


n  = 


u 

e'ToRo 

e‘YoRl 

4wiuie0 

Ro 

R1 

»0Cz+h) 


Jo(Xp)lT"dx 

u0 


(17) 


and 


n 


ll  cos  $  3  f  2(ul  ‘  uo)  -uQ (z+h)  t 

=  '  ^  *  3PJ0  Vfu0  ♦  Y§u7 


z  47riwe 


( Xp ) XdX  , 


(18) 
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where  R2  =  p2  +  (z  -  h)2,  R2  =  p2  *  (z  +  h)2,  u2  =  A2  +  y2,  and  u2  =  a2  +  y 2, 
From  equations  (17)  and  (18), 

J  *  H  cos  «  3  e'Y°R0  e'Y°Rl  f°2yle 

4™^  *  3p  Rq  -  Rl  *J  y2UQ  +  y2UiJ0^^X  .  (19 


3 

e"Y  0R0 

e'Y0Rl 

3p 

Ro 

R1 

-J0(Xp)AdXf.  (19) 


If  Jyqz1  |  £  0.5  (i.e.,  |  \Jn2  -  1 1  >_  2),  which  is  applicable  in  most  prac¬ 
tical  cases  (for  example,  |n2|  >_  81  at  all  frequencies  for  the  sea/air  case). 


ui  -  U0  -Und 


U1  +  uo 


/U!  '  Uo\  _  2U0 

’  \U1  +  uo/  "  U1  +  uc 


-u.d 

1  -  e  0 


«  1), 


For  small  Sommerfeld  numerical  distances  (i.e.,  |y  R  Vn2  -  l/(2n3) 


y5uo  +  Youi  '  Y?u0  • 


Utilization  of  the  identity  (Uj  -  u0)(Uj  +  uQ)  =  y2  -  y2,  equations  (21) 
and  (22),  and  Sommerfeld ' s  integrallO 


CO 

P  .[  e  “0 

Jo 


(z+h)T  „  ,  A,,  e'Y0Rl 

JqCAp^A  = 


results  in 


x  *  4iriwe0  I  RQ 


U  je  ~Y°  R°  e-YoRi  e-YoR:  e-YoR2 

R0  R1  R1  R2 


U  e 


4 7T  i. (x>£  o  I  Rq 


_Y0R0  'Y0R2 


■t  %  1 l  cos 

v  •  n  -  — : - -  x 

4iriue0 


I  ftp  cos 
4Trioie„ 


« |  \  "7  R>  I 

L|(1  * 


e’YoRi 

*  Y»R>)— 
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n 


In  COS  jl 
z  '  dTiwe. 


« -  (*  - 


in  cos  $ 
4Tricue„ 


f  /  -u.,d\  -u.(z+h)  l 

J  ^1  -  e  0  Je  01  "Jjfxp)- 


(1  -  1/n2)  /  ll  -  e 


C  26 ) 


2 

0 


where  R2  =  P2  +  (d  +  z  +  h) 2 . 


Because  of  the  u2  term  in  the  denominator  of  equation  (26) ,  it  does  not 
readily  appear  that  the  equation  can  be  expressed  in  closed  form,  except  for 
the  quasi-static  range  (y0  ~  0) ,  where  equation  (26)  reduces  to 


„  in  cos  ij )  f  (,  -Xd^  -X(z+h) 

n-  '  -  4nime„  J  V  *  6  /e  Jl(Xp)dX 


in  cos  <j>  |(d  +  z  +  h)  _  (z  +  h) 


47ri(jjeQp 


kl 


27) 


However,  if  we  replace  X2  by  u2  -  Yg.  equation  (26)  can  be  broken  up  nto 
four  integrals,  two  of  which  are  of  the  type 


00 

/  -unz  1 

-Ynz  Z  *YnR 

1  e  0  J1(Xp)dX  =  ~ 

6  K 

(28) 


where  R2  =  p2  +  z2,  and  two  of  which  are  of  the  type 


^0 


_0O 

j.  3  1 


(Xp)dX  . 


o  o 


(29) 


Since 


12 


a 

/ 


-uoz 


6  U-VXp)dA  =  77 


-  e‘Y°R 


(30) 


then 
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/-u  z 

e  1  j 
“2  1 


-Ynz 


(Ap)dA  = 


■n 


-YnR. 


dz  . 


(31) 


Inserting  equations  (28)  and  (31)  into  equation  (26)  results  in 


IS,  cos  $(1  -  1/n2) 

(d  +  z  +  h) 

4iriwe0p 

r2 

R1 


where 


(32) 


I  -  -Yr 


[  .  e’WRl> 

* 

dz 


When  (z  +  h)  »  p,  -  (z  +  h)  ,  and  R2  -  R  -  d,  resulting  in 

i  -  (i  -  e-T°d)e'T°(2*h)  . 


For  | yd |  <  0.5  (i.e.,  |n2|  >  IS), 


t  j  "Vz+h) 

I  '  Y0de  u 


(33) 


(34) 


(35) 


Comparing  equations  (27)  and  (32)  ,  we  see  that  for  the  quasi-static  range 
( I Y0rx I  <<:  1) »  I  is  negligible  compared  to  the  other  two  terms  in  equation 
(32).  Therefore,  if  we  assume  from  the  outset  that  Rj  >>  |d|,  equation  (33) 
becomes 


I  -  yQde 


*Y0Rl 


(36) 


which  is  identical  to  equation  (35)  when  (z  +  h)  >>  p. 
Therefore,  we  can  approximate  the  HED  n  vector  as 


IS,  cos  4> ( 1  -  1/n2) 

_  Y  R 

(d  +  z  +  h)e  02 

4iriu)e0p 

r2 

(z  +  h)e  Y°Rl  ,  *YnR j 

1 — i - +  Ve 


(37) 


subject  to  the  restriction  |n2|  >  15.  This  equation  should  also  be  valid  for 
|n2|  >  5  if  Ygd  is  replaced  by  (l  -  e  YqC*)  . 

Mow,  it  seems  that  we  have  gone  to  a  lot  of  trouble  to  derive  the  last 
term  of  the  n  expression.  However,  it  is  this  term  that  yields  the  vertically 
polarized  farfield  (for  small  numerical  distances)  produced  by  a  HED  located 
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very  near  or  below  the  earth's  surface.  This  is  the  so-called  quadripole  term 
(i.e.,  the  quadripole  moment  is  Iiy  d  =  2l£/n) .  Because  of  the  close  spacing 
and  opposite  sense  of  the  dipole  and  its  image,  direct  radiation  from  the 
dipole  is  not  the  prime  mechanism.  Rather,  the  quadripole  consisting  of  ver¬ 
tical  conduction  currents  in  the  lossy  medium  represents  the  prime  source. 13 


If  the  first  two  terms  in  equation  (37)  are  ignored  and  z  =  h  =  0,  equa¬ 
tion  (37)  reduces  to 


n 


z 


Ii.  COS  (j>  (  ,  YnP 

-t~. - r|Y„de  u 

4iriu)£0P  \'0  , 


(38) 


When  |  yqp  !  >>  1, 


ito£ 


an 

_ z 

o  ap 


H  COS  0,  2\ 
2nr1p  VO/ 


“V 


(39) 


which  is  the  correct  farfield  result  for  small  numerical  distances  when  |n2| 
»  1. 


Since  we  have  now  derived  expressions  for  the  HED  Hertz  vector  (equations 
(24),  (25  .  and  (37)),  the  fields  in  air  can  be  obtained  from 


?  =  -y^n  +  $($  •  n) 

H  =  iue0(V  x  jf)  , 


(40) 


The  ssulting  HED  finitely  conducting  earth-image  theory  field  expressions 
for  the  ai  -to-air  propagation  case  are  presented  in  table  1.  They  are  valid 
for  small  numerical  distances  and  |n2j  >  15.  When  |y0riI  «  1 ,  they  reduce  to 

the  quasi-static  range  image-theory  results. When  |y0riI  >?  1*  they  reduce 
to  Norton's^.  15  farfield  results. 


These  results  easily  can  be  extended  to  a  multilayered  earth  simply  by 
letting  d  =  (2/y1)Q»  where  Q  is  the  familiar  plane-wave  correction  factor 
employed  to  account  for  the  presence  of  stratification  in  the  earth. 1&.17 


HMD  AIR-TO-AIR  PROPAGATION 


Four  of  the  six  HMD  expressions  valid  for  h  and  z  >_  0  can  be  obtained 
completely  from  reciprocity  consideration  (Ep,  E^,  Ez,  and  H2) .  The  remaining 
two  can  be  obtained  from  Maxwell's  equations  (Hp  and  H^) .  Alternatively,  they 
can  be  obtained  from 
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Table  1.  HED  Air-to-Air  Propagation  Equations. 


■T 


TR  6511 


U  cos  <t> 

/3p2 

■V  *  v0R„) 

-YoRo 

t  ..2-2  e 

'i  -4 

1  • 

1  •. 

o  4iriajeg 

h ' 

°  J  R^ 

1 

•  (1  +  YgRj)  +  YgP2 

-Y0Rl 

e  ..2 

'e~YoR° 

e"Y0R2 

R?  *0 

1 

Ro 

r2  j 

10 


'YoR0 


4Tri(jL)£, 


o  *  (l  -  4r)»  *  r„R,^ 

0 


‘YqR1 


+  Y 


-YOr0  -Y0R2 
2|e _ e _ 


0\  R, 


)l 


Up  cos  <{> 

Cz  -  h)  / 

4iriojeQ 

Rs  V 
0 

(3  ♦  3yqR0  +  y2R2)e 


-  (>  -  ^■)l£^j!i(3  *  *.*,  * 


YoR0 
-Y0Ri 


(d  +  2  +  h)  -YqR2  +  Y0°  I 


YndR,  -  (z  +  h) 


-Y0r1 


U  sin  4> 

(d  ♦  z  +  h) 

4ir 

R2 

Rl 


-y0R2 


YgdRj  -  (z  +  h) 


p2Ri 


-y0r!  .  +  YoVe-YoR° 


*1 


U  cos  4> 

(Z  -  h). 

4tt 

D  3  1 

Ro 

Rl 


Cd  +  z  +  h) _'Y0R2 

(z  +  h)  -y„Ri 

.  de'Y°Ri 

V  R  + 

R2 

R1  J 

+  P2Rj 

YoKi  Y0P 

,,  Up  sin  <j> 

n  - 1 — - - 

z  4tt 


-Y0Ro 


(1  ♦  YnRn) 


0  0'  d 3 
K0 


-  C1  +  Y„R,) 


'y0R2 


0  2'  R | 


■V 


TR  6511 


Ho  =  'Yony  sin  *  + 


H4>  =  -YSny  cos  ♦  +  j  x  ^  *  “)  • 


Following  the  same  procedure  outlined  in  the  derivation  of  the  HED  field 
components  results  in 


IA 

e"Y°R° 

4tt 

Ro 

-y„R„ 

*  *  IAp  sin  d>  e  0 

*  •  11  '  ‘  •  1  V- '  r3  1  +  YoV  + 

0 


-a  +  y„r,j 


The  resulting  HMD  finitely  conducting  earth-image  theory  field  expressions 
for  the  air-to-air  propagation  case  are  presented  in  table  2.  They  are  valid 
for  | n2 |  >  15  and  small  numerical  distances.  When  Iy0riI  <<  1,  they  reduce  to 
the  quasi-static  range  image-theory  results. *»2-  When  |y0Rj|  >>  they  reduce 
to  Norton's14 farfield  results. 


SUBSURFACE-TO-AIR  PROPAGATION 


The  HED  and  HMD  image-theory  expressions  for  the  subsurface-to-air  propa¬ 
gation  case  (h  <_  0,  z  ^  0)  can  be  obtained  from  the  air-to-air  propagation 
equations  (tables  1  and  2)  simply  by  setting  h  =  0  and  multiplying  each  expres¬ 
sion  by  exp(Yjh).  The  resulting  equations  are  presented  in  tables  3  and  4  and 

will  be  valid18  for  R  =  ^p^  +  z2-  >  |3hj,  |n2|  >  15,  and  small  numerical  dis¬ 
tances.  It  should  be  noted  that  by  following  the  procedure  outlined  by  Bannis¬ 
ter  and  Dube,18  the  restriction  R  >  |3h|  can  become  less  stringent. 

When  JYiR|  >>  1,  the  HED  and  HMD  subsurface-to-air  propagation  equations 

reduce  to  the  nearfield  and  farfield  range  results  presented  in  tables  3.1  and 
3.3  of  Kraichman19  (see  also  Bannister2). 

AIR-TO-SUBSURFACE  PROPAGATION 


The  HED  and  HMD  image-theory  expressions  for  the  air-to-subsurface  propa¬ 
gation  case  (h  >_  0,  z  <_  0)  can  be  obtained  from  the  air-to-air  propagation 
equations  (tables  )  and  2)  simply  by  setting  z  ■  0  and  multiplying  each  expres¬ 
sion  by  exp(YjZ).  (Both  Ez  components  must  also  be  multiplied  by  1/n2  to 
satisfy  the  boundary  conditions.)  The  resulting  equations  are  presented  in 
tables  5  and  6  and  will  be  valid18  for  R'  *  ^p2  +  h2  >  |3z|,  |n2|  >  15,  and 
small  numerical  distances. 
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Table  2.  HMD  Air-to-Air  Propagation  Equations 
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Table  3  HED  Subsurface-to-Air  Propagation  Equations 
[r2  =  p2  +  z2,  r2  =  p2  +  (d  +  Z)2]. 
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Table  4.  HMD  Subsurface-to-Air  Propagation  Equations 
[R2  =  p2  +  z2,  R2  =  p2  +  (d  +  z)2] 
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Table  5.  HED  Air-to-Subsurface  Propagation  Equations 
[(R')2  =  p2  ♦  h2,  (R’)2  =  p2  +  (d  +  h)2] 
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Table  6.  HMD  Air-to-Subsurface  Propagation  Equations 
[(R')2  =  p2  +  h2,  (R? ) 2  =  p2  +  (d  +  h)2] 
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When  |yjR' |  >>  1,  the  HED  and  HMD  air-to-subsurface  propagation  equations 
reduce  to  the  nearfield  range  results  presented  in  table  3.5  of  Kraichman. 


SUBSURFACE-TO-SUBSURFACE  PROPAGATION 


The  HED  and  HMD  image-theory  expressions  for  the  subsurface-to-subsurface 
propagation  case  (h  <_  0,  z  <_  0)  can  be  obtained  from  the  air-to-air  propagation 
equations  (tables  1  and  2)  simply  by  setting  both  z  and  h  equal  to  zero  and 
multiplying  each  expression  by  exp^Cz  +  h)  ] .  (Both  Ez  components  must  also 
be  multiplied  by  1/n2  =  y2/y2  to  satisfy  the  boundary  conditions.)  The 
resulting  expressions  are  presented  in  tables  7  and  8  and  should  be  valid18 
(for  most  cases)  for  p  >  3 | z  +  h|,  |n2|  >  15,  and  small  numerical  distances. 
When  |y0p|  <<  1  and  z  =  h  =  0,  these  expressions  reduce  to  the  quasi-static 
range  surface-to-surface  image-theory  results. 1*2  When  !y,p|  >>  1,  they  reduce 
to  the  subsurface-to-subsurface  nearfield  and  farfield  range  results  presented 
in  tables  3.2  and  3.7  of  Kraichman. 


COMPARISON  WITH  EXACT  SOMMERFELD  INTEGRATION  RESULTS 


Mittra  et  al.20  have  presented  some  exact  Sommerfeld  integration  results 
for  the  HED  0nx  and  ll2  vectors  (i.e.,  the  correction  terms  to  the  perfectly 
conducting  ground  solution)  for  frequencies  of  3  to  30  MHz.  For  their  case, 

Rj  =  10  m,  9  =  tan_1[p/(z  +  h)  ]  =  10°,  <p  =  0°,  and  the  quantity  IQ  =  If/(iu»e0) 
is  normalized  to  unity. 


The  image-theory  solution  of  n2  is  given  by  equation  (37) ,  while  QIIX  is 
given  by  the  last  two  terms  of  equation  (24);  that  is. 


0 


4ir 


-YoR 


0  1 


(44) 


A  comparison  of  the  image  theory  and  exact  Sommerfeld  integration  results 
for  the  situation  where  er  =  40  and  Oj  =  1  S/m  is  presented  in  figure  3.  For 

this  situation,  n2  varies  from  40  -  j 150  at  3  MHz  to  40  -  j 15  at  30  MHz.  Since 
|n2|  >  15,  d_.  -  d__  ~  d  =  2/y , .  Note  that  the  agreement  between  the  two  solu- 

1 M  It  i 

tions  is  excellent. 


Presented  in  figure  4  is  a  comparison  of  the  image-theory  and  exact  Som¬ 
merfeld  integration  results  for  the  case  where  er  =  10  and  Pj  =  10-2  S/m.  For 

this  case,  n2  varies  from  10  -  j6  at  3  MHz  to  10  -  j0.6  at  30  MHz.  Since 
n2 :  <  15  and  (z  ♦  h)  >>  p,  we  have  replaced  Ygd  by  1  -  exp(-YQd)  in  the  image 
theory  JI  expression  (equation  (37)).  Note  that  the  agreement  is  excellent 
(within  1  percent  for  the  0nx  component  and  very  good  (within  5  percent)  for 
the  component. 
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Table  7.  HED  Subsurface-to-Subsurface  Propagation 
Equations  (p?  =  p2  +  d2) 
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Table  8.  HMD  Subsurface-to-Subsurface  Propagation 
Equations  (p2  =  p2  +  d2) 
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EXTENSION  TO  LARGE  NUMERICAL  DISTANCES 


Image  theory  also  can  be  utilized  to  determine  the  fields  at  large  numer¬ 
ical  distances  (i.e.,  |p|  =  |-YQR1/2n2j  >>  1).  For  the  sake  of  simplicity,  we 
will  let  z  =  h  =  0  and  |n2|  >>  1.  For  this  case,  equation  (18)  reduces  to 
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where  D2  =  p2  +  (n2d/2)2. 

Because  p2  »  (n2d/2) 2 , 
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which  is  the  correct  farfield  result  for  large  numerical  distances. 
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CONCLUSIONS 


Simple  engineering  expressions  for  HED  and  HMD  air-to-air,  subsurface-to- 
air,  air-to-subsurface ,  and  subsurface-to-subsurface  propagation  have  been 
derived  by  employing  finitely  conducting  earth-image  theory  techniques.  For 
the  air-to-air  propagation  case,  the  expressions  are  valid  from  the  quasi-static 
to  the  farfield  ranges  as  long  as  |n2l  >  15  and  the  Sommerfeld  numerical  dis¬ 
tance  is  small.  For  the  subsurface-to-air ,  air-to-subsurface,  and  subsurface- 
to-subsurface  cases,  the  additional  restriction  that  the  measurement  distance 
be  greater  than  three  times  the  burial  depth  of  the  source  and/or  receiver 
must  be  met.  We  have  also  demonstrated  that  image  theory  can  be  utilized  to 
determine  the  fields  at  large  numerical  distances. 

We  have  compared  successfully  image-theory  and  exact  Sommerfeld  integra¬ 
tion  results  for  four  cases,  yielding  agreement  within  1  percent  for  three 
comparisons  and  within  5  percent  for  the  other. 

It  should  be  noted  that  the  two  media  can  be  inverted  and  the  air 
replaced  by  the  earth's  crust  (of  conductivity  a2  and  dielectric  constant  e2) . 
The  same  equations  (tables  1  through  8)  can  be  utilized  as  long  as 
|y2R1/(2n2)|  <<  1  and  | n| |  =  lyf/yfl  >  *5  simply  by  replacing  iwe0  by  a2  +  iu)e2 . 

The  results  presented  in  this  report  should  be  particularly  useful  for 
sea/air  and  sea/earth' s -crust  propagation.  They  should  also  be  helpful  to 
geophysicists  engaged  in  determining  the  electrical  properties  of  the  earth: 

The  simple,  yet  accurate,  formulas  obtained  from  this  theory  make  it  a  very 
strong  tool  and  a  promising  one  for  determining  the  coupling  between  antennas 
located  above  or  below  the  earth's  surface. 
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